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Abstract 

Helical configurations of inhomogeneous symmetric rods with non-constant bend- 
ing and twisting stiffness are studied within the framework of the KirchhofF rod 
model. Prom the static Kirchhoff equations, we obtain a set of differential equations 
for the curvature and torsion of the centerline of the rod and the Lancret 's theorem 
is used to find helical solutions. We obtain a free standing helical solution for an 
inhomogeneous rod whose curvature and torsion depend on the form of variation of 
the bending coefficient along the rod. These results are obtained for inhomogeneous 
rods without intrinsic curvature, and for a particular case of intrinsic curvature. 
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1 Introduction 



Helical filaments are tridimensional structures commonly found in Nature. 
They can be seen in microscopic systems, as biomolecules [1], bacterial fibers [2] 
and nanosprings [3], and in macroscopic ones, as ropes, strings and climbing 
plants [4,5,6]. Usually, the axis of all these objects is modeled as a circular 
helix, i. e. a 3D-space curve whose mathematical geometric properties, namely 
the curvature, kp, and the torsion, r F , are constant [7,8,9]. This kind of helical 
structure has been shown to be a static solution of the Kirchhoff rod model [7]. 

The Kirchhoff rod model [10,11] has been proved to be a good framework to 
study the statics [7,12,13] and dynamics [29] of long, thin and inextensible 
elastic rods. Applications of the Kirchhoff model range from Biology [1,14,5] 
to Engineering [15] and, recently, to Nanoscience [16]. In most cases, the rod or 
filament is considered as being homogeneous, but the case of nonhomogeneous 
rods have also been considered in the literature. It has been shown that non- 
homogeneous Kirchhoff rods may present spatial chaos [17,18]. In the case of 
planar rods, Domokos and collaborators have provided some rigorous results 
for non-uniform elasticae [19] and for constrained Euler buckling [20,21]. Devi- 
ations of the helical structure of rods due to periodic variation of the Young's 
modulus were verified numerically by da Fonseca, Malta and de Aguiar [22]. 
Nonhomogeneous rods subject to given boundary conditions were studied by 
da Fonseca and de Aguiar in [23]. The effects of a nonhomogeneous mass 
distribution in the dynamics of unstable closed rods have been analyzed by 
Fonseca and de Aguiar [24]. Goriely and McMillen [25] studied the dynamics 
of cracking whips [26] and Kashimoto and Shiraishi [27] studied twisting waves 
in inhomogeneous rods. 
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The stability analysis of helical structures is of great importance in the study 
of the elastic behavior of filamentary systems and has been performed both 
experimentally [28] and theoretically [29,30,31]. It has been also shown that 
the type of instability in twisted rods strongly depends on the anisotropy of 
the cross section [32,33,34]. 

Here, we consider a rod with nonhomogeneous bending and twisting coef- 
ficients varying along its arclength s, B(s) and C(s), respectively. We are 
concerned with the following question: is there any helical solution for the sta- 
tionary Kirchhoff equations in the case of an inhomogeneous rod ? The answer 
is 'yes' and it will be shown that the helical solution for an inhomogeneous 
rod with varying bending coefficient cannot be the well known circular helix, 
for which the curvature, hp, and torsion, t>, are constant. To this purpose, 
we shall derive a set of differential equations for the curvature and the torsion 
of the centerline of an inhomogeneous rod and then apply the condition that 
a space curve must satisfy to be helical: the Lancret's theorem. We shall ob- 
tain the simplest helical solutions satisfying the Lancret's theorem and show 
that they are free standing helices, i.e., helices that are not subjected to axial 
forces [30] . A resulting helical structure different from the circular helix, from 
now on, will be called a Lancret helix. 

According to the fundamental theorem for space curves [9], the curvature 
kp(s), and the torsion, t f (s), completely determine a space curve, but for 
its position in space. We shall show that the &f(s) and Tf(s) of a Lancret he- 
lix depend directly on the bending coefficient, B(s), an expected result since 
the centerline of the rod does not depend on the twisting coefficient (see for 
example, Neukirch and Henderson [13]). 
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Some motivations for this work are related to defects [35] and distortions [36] 
in biological molecules. These defects and distortions could be modeled as 
inhomogeneities along a continuous elastic rod. 



In Sec. II we review the general definition of a space curve, the Frenet basis 
and the so-called Lancret's theorem. In Sec. Ill we present the static Kirchhoff 
equations for an intrinsically straight rod with varying stiffness, and derive 
the differential equations for the curvature and torsion of the rod. In Sec. 
IV we use the Lancret's theorem for obtaining helical solutions of the static 
Kirchhoff equations and we show that they cannot be circular helices if the 
bending coefficient is not constant. As illustration, we compare a homogeneous 
rod with two simple cases of inhomogeneous rods: (i) linear and (ii) periodic 
bending coefficient varying along the rod. The circular helix has a well known 
relation of the curvature and torsion with the radius and pitch of the helix. In 
Sec. V we define a function involving all these variables in such a way that for 
the circular helix its value is identically null. We have verified, numerically, 
that this function approaches zero for the inhomogeneous cases considered 
here. In Sec. VI we analyse the cases of null torsion (straight and planar 
rods). Since helical solutions of intrinsically straight rods are not dynamically 
stable [30], in Sec. VII we consider a rod with a given helical intrinsic curvature 
and we obtain, for this helical solution of the static Kirchhoff equations 

similar to that of an intrinsically straight inhomogeneous rod. In Sec. VIII we 
summarize the main results. 
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2 Curves in space 

A curve in space can be considered as a path of a particle in motion. The 
rectangular coordinates (x, y, z) of the point on a curve can be expressed as 
function of a parameter u inside a given interval: 

x = x(u) , y = y{u) , z = z(u) , u\ < u < u<i . (1) 

We define the vector x(u) = (x(u),y(u), z(u)). If u is the time, x(tt) represents 
the trajectory of a particle. 

2.1 The Frenet frame and the Frenet-Serret equations 

The vector tangent to the space curve at a given point P is simply dx./du. It 
is possible to show [9] that if the arclength s of the space curve is considered 
as its parameter, the tangent vector at a given point P of the curve x(s) is 
a unitary vector. So, using the arclength s to parametrize the curve, we shall 
denote by t its tangent vector 

dx. , , 

||t|| = 1. The tangent vector t points in the direction of increasing s. 

The plane defined by the points Pi, P 2 and P 3 on the curve, with P 2 and P 3 
approaching P 1 , is called the osculating plane of the curve at Pi [9]. Given a 
point P on the curve, the principal normal at P is the line, in the osculating 
plane at P, that is perpendicular to the tangent vector at P. The normal 
vector n is the unit vector associated to the principal normal (its sense may 
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be chosen arbitrarily, provided it is continuous along the space curve). 

From t.t = 1, differentiating with respect to s (indicated by a prime) it follows 
that: 



t.t' = , 



(3) 



so that t and t' are orthogonal. It is possible to show that t' lies in the 
osculating plane, consequently t' is in the direction of n. This allows us to 
write 



hp being called the curvature of the space curve at P. 

The curvature measures the rate of change of the tangent vector when moving 
along the curve. In order to measure the rate of change of the osculating plane, 
we introduce the vector normal to this plane at P: the binomial unit vector 
b. At a point P on the curve, b is defined in such a way that 



The frame {n, b, t} can be taken as a new frame of reference and forms the 
moving trihedron of the curve. It is commonly called the Frenet frame. The rate 
of change of the osculating plane is expressed by the vector b'. It is possible 
to show that b' is anti-parallel to the unit vector n [9]. So we can write 



t' = kF n , 



(4) 



b = t x n . 



(5) 



b' = - 



(6) 



tf being called the torsion of the space curve at P. 



6 



The rate of variation of n [9] can be obtained straightforwardly. It is given by 
n' = —kp t + rpb . (7) 



The set of differential equations for {t, n, b} is 



t' = kF n , 



n' = —k F t + r F b , 



(8) 



b' = — t f n , 

and are known as the formulas of Frenet or the Serret-Frenet equations [9]. 
2.2 The Fundamental theorem of space curves 

A space curve parametrized by its arclength s is defined by a vectorial func- 
tion x(s). The form of x(s) depends on the choice of the coordinate system. 
Nevertheless, there exists a form of characterization of a space curve given 
by a relation that is independent of the coordinates. This relation gives the 
natural equation for the curve. 

/cf(s) gives the natural equation in the case of planar curves. Indeed, if ip is 
the angle between the tangent vector of the planar curve and the rr-axis of 
the coordinate system, it is possible to show that k,F = dtp/ds. Since cos(^) = 
dx/ds and sin((^) = dy/ds, knowing then (p(s), x(s), and y(s) of the 

planar curve can be obtained immediatly: 



s s 



s 



ip(s) = / kp(s)ds, x(s) = / cos ip(s)ds, 
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In the case of non-planar curves, if we have two single valued continuous func- 
tions fci?(s) and tf(s), s > 0, then there exists one and only one space curve, 
determined but for its position in space, for which s is the arclength, kp(s) 
the curvature, and tf(s) the torsion. It is the Fundamental theorem for space 
curves [9]. The functions kp(s) and r F (s) provide the natural equations of the 
space curve. 



2.3 Curves of constant slope: the Lancret's theorem 

A space curve x(s) is a helix if the lines tangent to x make a constant angle 
with a fixed direction in space (the helical axis) [8,9]. Denoting by a the unit 
vector of this direction, a helix satisfies 

t.a = cos a = constant . (10) 

Differentiating Eq. (10) with respect to s gives a.n = 0. Therefore a lies in 
the plane determined by the vectors t and b: 

a = t cos a + b sin a . (11) 

Differentiating Eq. (11) with respect to s, gives 
= (hp cos a — t f sina)n , 



or 



kp 

— = tana = constant. (12) 



8 



This result says that for curves of constant slope the ratio of curvature over 
torsion is constant. Conversely, given a regular curve for which the equation 
(12) is satisfied, it is possible to find [9] a constant angle a such that 

n (k F cos a — t f sin a) = , 

d , . 
— (t cos a + b sin a) = , 
as 

implying that the vector a = t cos a + b sin a is the unit vector along the axis. 
Moreover, a.t = cos a=constant, so that the curve has constant slope. This 
result can be expressed as: 

A necessary and suficient condition for a space curve to be a curve of constant 
slope (a helix) is that the ratio of curvature over torsion be constant. It is the 
well known Lancret's theorem, dated of 1802 and first proved by B. de Saint 
Venant [9,37]. 

If a helical curve x(s) is projected onto the plane perpendicular to a, the 
vector xi(s) representing this projection is given by 



It is possible to show [9] that the curvature k\ of the projected curve is given 



The shape of the planar curve obtained by projecting a helical curve onto the 
plane perpendicular to its axis is used to characterize it. For example, the well 
known circular helix projects a circle onto the plane perpendicular to its axis. 




(13) 



by: 




(14) 
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The spherical helix projects an arc of an epicycloid onto a plane perpendicular 
to its axis [9]. The logarithmic spiral is the projection of a helical curve called 
conical helix [9]. 



3 The static Kirchhoff equations 

The statics and dynamics of long and thin elastic rods are governed by the 
Kirchhoff rod model. In this model, the rod is divided in segments of infinites- 
imal thickness to which the Newton's second law for the linear and angular 
momentum are applied. We derive a set of partial differential equations for 
the averaged forces and torques on each cross section and for a triad of vectors 
describing the shape of the rod. The set of PDE are completed with a linear 
constitutive relation between torque and twist. 

The central axis of the rod, hereafter called centerline, is represented by a 
space curve x parametrized by the arclength s. A Frenet frame is defined for 
this space curve as described in the previous section. For a physical filament 
the use of a local basis, {d 1; d 2 , d 3 }, to describe the rod has the advantage of 
taking into account the twist deformation of the filament. This local basis is 
defined such that d 3 is the vector tangent to the centerline of the rod (d 3 = t), 
and d x and d 2 lie on the cross section plane. The local basis is related to the 
Frenet frame {n, b,t} through 



\ 



cos £ — sin £ 



(di d 2 dg) = (n b t) 



sin £ cos £ 



(15) 



V 



o 







1 
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where the angle £ is the amount of twisting of the local basis with respect to 
t. 



In this paper, we are concerned with equilibrium solutions of the Kirchhoff 
model, so our study departs from the static Kirchhoff equations [38] . In scaled 
variables, for intrinsically straight isotropic rods, these equations are: 

F' = , (16) 
M' = Fxd 3 , (17) 
M = B(s) h dx + B(s) k 2 d 2 + C(s) k 3 d 3 , (18) 

the vectors F and M being the resultant force, and corresponding moment 
with respect to the centerline of the rod, respectively, at a given cross section. 
As in the previous section, s is the arclength of the rod and the prime ' denotes 
differentiation with respect to s. hi are the components of the twist vector, 
k, that controls the variations of the director basis along the rod through the 
relation 

d| = k x di , i = 1,2,3 . (19) 



k\ and k 2 are related to the curvature of the centerline of the rod (/jp = 
\Jk\ + fcf) and k 3 is the twist density. B(s) and C(s) are the bending and 
twisting coefficients of the rod, respectively. In the case of macroscopic fila- 
ments the bending and twisting coefficients can be related to the cross section 
radius and the Young's and shear moduli of the rod. Writing the force F in 
the director basis, 

F = /idi + f 2 d 2 + / 3 d 3 , (20) 



the equations (16-18) give the following differential equations for the compo- 
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nents of the force and twist vector: 



f[-f2h + hk 2 = 0, (21) 
f 2 + fik 3 -f 3 h = 0, (22) 

/3-/ifc 2 + / 2 A;i = 0, (23) 

(B(s) h)> + (C(s) - B(s)) k 2 k 3 -f 2 = 0, (24) 

(B(s) k 2 )' - (C(s) - B(s)) h k 3 + h = , (25) 

(C(s) k 3 )' = . (26) 

The equation (26) shows that the component M 3 = C(s)k 3 of the moment 
in the director basis (also called torsional moment), is constant along the 
rod, consequently the twist density k 3 is inversely proportional to the twisting 
coefficient C(s) 

M«) - ^ • (27) 

In order to look for helical solutions of the Eqs. (21-26) the components of 
the twist vector k are expressed as follows: 

k 1 = k F (s)sin^, (28) 
k 2 = k F (s) cos£ , (29) 
k 3 = e + r F (s), (30) 



where k F (s) and t f (s) are the curvature and torsion, respectively, of the space 
curve that defines the centerline of the rod and £ is given by Eq. (15). If the 
rod is homogeneous, the helical solution has constant k F and t f , and £' is 
proved to be null [5]. 

Substituting Eqs. (28-30) in Eqs. (21-26), extracting f\ and f 2 from Eqs. (25) 
and (24), respectively, differentiating them with respect to s, and substituting 
in Eqs. (21), (22) and (23), gives the following set of nonlinear differential 
equations: 
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[M 3 k F (s) " B ( s ) M*) r F (s)]' - (B(s) k F (s))' r F (s) = , (31) 
(B(s) k F (s))" + k F (s) r F (s)[M 3 - B(s) r F (s)} - f 3 (s) k F (s) = , (32) 

(B(s)k F (s))'k F (s) + ft(s) = 0. (33) 



Appendix A presents the details of the derivation of Eqs. (31-33). 

The Eqs. (31-33) for the curvature, k F , and torsion, r F do not depend on the 
twisting coefficient, C(s). Therefore, the centerline of an inhomogeneous rod 
does not depend on the twisting coefficient like in the case of homogeneous 
rods (see, for example, Eqs. (13) and (14) of Ref. [13]). 

Langer and Singer [39] have obtained a set of first-order ordinary differential 
equations for the curvature and torsion of the centerline of a homogeneous 
rod that contains terms proportional to k\ and rp. The Eqs. (31-33) have the 
advantage of involving only terms linear in k F and r F . 

4 Helical solutions of inhomogeneous rods 

In order to find helical solutions for the static Kirchhoff equations, we apply 
the Lancret's theorem to the general equations (31-33). We first rewrite the 
Lancret's theorem in the form: 



Ms) = pT F (s) , 



(34) 



with p^O. From Eq. (12), 



(3 = tan a = Constant . 



(35) 
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Substituting Eq. (34) in Eq. (31) we obtain 

t' f (M 3 -Bt f )-2t f (B t f )' = . (36) 

Substituting Eq. (34) in Eq. (32) and extracting / 3 , we obtain 

U = + r F (M 3 -Bt f ). (37) 

Differentiating / 3 with respect to s and substituting in Eq. (33) we obtain the 
following differential equation for t f : 



{Bt f )"> (Bt f )"t 



F + ((3 2 + l)T F (Br F y = 0, (38) 



t f t 2 f 



where the Eq. (36) was used to simplify the above equation. One immediate 
solution for this differential equation is 

(Bt f )' = 0, (39) 

that substituted in Eq. (36) gives 

t' f (M 3 - Bt f ) = . (40) 

For non-constant t F) the Eq. (40) gives the following solution for t f : 

T ' (s) = w, ■ < 41 » 

Substituting the Eqs. (39) and (41) in Eq. (37) we obtain that 

fs(s) = . (42) 
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Substituting Eq. (41) in (34) we obtain: 



k F (s) = (3 



M 3 



(43) 



B(s) ■ 



Substituting Eq. (15) in Eq. (20), the force F becomes 



F = (/xcosf -/ 2 sinf)n+ (/isinf + / 2 cosOb + / 3 t , 



(44) 



where {n, b, t} is the Frenet basis. Using the Eqs. (A. 7) and (A. 8) for fi and 
f'2 (Appendix A), we obtain 



where / 3 , in the inhomogeneous case, must satisfy the Eq. (33). 

Substituting the Eqs. (41-43) in the Eq. (45), and using Eq. (35), it follows 
F = 0. Therefore, the helical solutions satisfying (39) are free standing. 

Now, we prove that a circular helix cannot be a solution of the static Kirchhoff 
equations for a rod with varying bending stiffness. If a helix is circular, k' F = 
and t' f = 0, and from Eq. (31) we obtain: 



Since B'(s) 7^ 0, Eq. (46) will be satisfied only ii k F = and/or r F = 0. 
Therefore, it is not possible to have a circular helix as a solution for a rod 
with varying bending coefficient. 

The solutions for the curvature &f, Eq. (43), and the torsion tf, Eq. (41), can 
be used to obtain the unit vectors of the Frenet frame (by integration of the 
Eqs. (8)). From Eqs. (43), (35) and (11), we can obtain a and a once k F (0), 



F 



(Bk F )'n + k F [M 3 - B r F }b + f 3 t , 



(45) 



2k F r F B' = . 



(46) 
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M 3 and B(s) are given. By choosing the z-direction of the fixed cartesian basis 
as the direction of the unit vector a, we can integrate t in order to obtain the 
three-dimensional configuration of the centerline of the rod. 

Figure 1 displays the helical solution of the static Kirchhoff equations for rods 
with bending coefficients given by 



The case of constant bending (47) produces the well known circular helix 
displayed in Fig. la. Figs, lb-lc show that non-constant bending coefficients 
(Eqs. (48-49)) do not produce a circular helix. 

The helical solutions displayed in Fig. 1 satisfy the Lancret's theorem, Eq. (12). 
The tridimensional helical configurations displayed in Fig. 1 were obtained by 
integrating the Frenet-Serret equations (8) using the following initial condi- 
tions for the Frenet frame: t(s = 0) = (0, sin a, cos a), n(s = 0) = (—1,0,0) 
and b(s = 0) = (0, — cos a, sin a). This choice ensures that the z-axis is par- 
allel to the direction of the helical axis, vector a. The centerline of the helical 
rod is a space curve x(s) = (x(s),y(s), z(s)) that is obtained by integration 
of the tangent vector t(s). We have taken the helical axis as the z-axis and 
placed the initial position of the rod at x(0) = l/fci(0), y(0) = and z(0) = 
(in scaled units), where fci(0) is the curvature of the planar curve at s = 
obtained by projecting the space curve onto the plane perpendicular to the 
helical axis (Eq. (14)). From Eq. (14) we have 



Fig la: B a (s) = l, 

Fig lb: B b (s) = l + 0.007 s , 

Fig lc: £ c (s) = 1 + 0.1 sin(0.04s + 2) 



(47) 
(48) 
(49) 



fci(0) 




(50) 
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Using the Eq. (35), it follows that 

sin 2 a = — . (51) 



From Eq. (43), setting s = 0, we get 

p = tmm. ,52) 

M 3 

Substituting Eqs. (51) and (52) in Eq. (50), we obtain: 

x(0) = — = ^(°) g2 (°) f53 ) 



/cf(0) and M 3 are free parameters that have been chosen so that the helical 
solutions displayed in Fig. 1 have the same angle a. The parameters kp(0) = 
0.24 and M 3 = 0.05 give x(0) ~ 4 for the helical solutions displayed in Figs, 
la and lb, and the parameters &f(0) = 0.22 and M 3 = 0.05 give x(0) ~ 4.36 
for the helical solution displayed in the Fig. lc. 

For short the projection of the space curve onto the plane perpendicular to the 
helical axis will be called projected curve. As mentioned in Sec. II, the circle 
is the projected curve of the most common type of helix, the circular helix. 

Fig. 2 displays the projected curves related to the helical solutions displayed in 
Fig. 1. Fig. 2a shows that the helical solution of the inhomogeneous rod with 
constant bending coefficient projects a circle onto the plane perpendicular to 
the helical axis. 

If required, the natural equations for the projected curves displayed in Fig. 2 
are easily obtained, for instance, by substitution of the solution (Eq. (43)) for 
the curvature fcjr(s) of the helical rod into Eq. (14). The natural equation of 



17 



the projected curve is given by its curvature, 




(54) 



where (5 and sin" 2 a can be obtained by Eqs. (51) and (52). Then, in the 
Eq.(54), setting B(s) = -Bj(s), i — a,b, c, as given in Eqs. (47-49), produces the 
natural equation for the corresponding projected curve displayed in Fig. 2. The 
helical rod displayed in Fig. lb is a conical helix since the radius of curvature 
of its projected curve (inverse of fci(s)) is a logarithmic spiral (l/ki(s) is a 
linear function of s [9]). 

From Eqs. (28-30), (27) and (41) we obtain the variation of the angle £ between 
the local basis, d;, % — 1, 2, 3, and the Frenet frame, {n, b, t}: 



Eq. (55) shows that £' ^ for the general case of B(s) ^ C(s), i. e. helical 
filaments corresponding to inhomogeneous rods are not twistless. The circular 
helix is a helicoidal solution for the centerline of an inhomogeneous rod having 
constant bending coefficient. We emphasize that the inhomogeneous rod is not 
twistless in contrast with the homogeneous case where it has been proved that 



A homogeneous rod has B(s) and C(s) constant so that k 3 = Constant (from 
Eq. (26)). Since £' has been proved to be null for a helical solution of a ho- 
mogeneous rod (see reference [5]), Eq. (30) shows that the torsion t f must 
be a constant. In order to satisfy the Lancret's theorem (Eq. (12)) the curva- 
ture hp of the helical solution must also be a constant. Therefore, the only 
type of helical solution for a homogeneous rod is the circular helix, while an 




(55) 



r = 0[5]. 
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inhomogeneous rod may present other types of helical structures. 
5 Radius and Pitch of the helical solution 

The radius 1Z of a helix is defined as being the distance of the space curve to 
its axis. The pitch V of a helix is defined as the height of one helical turn, i.e., 
the distance along the helical axis of the initial and final points of one helical 
turn. 

For a circular helix, 71, V, k F and r F are constant, and it is easy to prove that 



For other types of helix, it constitutes a very hard problem in differential 
geometry to obtain the relation between the curvature k F and the torsion r F 
with the radius 7Z and the pitch V . We have seen in Sec. II that the definitions 
of curvature and torsion involve the calculation of the modulus of the tangent 
and normal vectors derivative with respect to the arclength of the rod. We also 
saw that the Frenet-Serret differential equations for the Frenet frame depend 
on the curvature and torsion. The difficulty of integration of the Frenet-Serret 
equations for the general case where k F (s) and r F (s) are general functions of 
s poses the problem of finding an analytical solution for the centerline of the 
rod, thus the difficulty of relating non constant curvature and torsion with non 
constant radius and pitch. Due to this difficulty we shall test the possibility 
of generalizing the relation (56) to the present inhomogeneous case. In order 
to do so, from the equation (56) we define: 




i 



(56) 




(57) 
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where TZ(s) and V(s) are the radius and the pitch of the helical structure as 
function of s. In the case of a circular helix, from Eq. (56), g\(s) = for all s. 

Since the z-axis is defined as being the axis of the helical solution we can 
calculate the radius TZ(s) through: 1Z(s) = \Jx 2 (s) + y 2 (s), where x(s) and 
y(s) are the x and y components of the vector position of the centerline of the 
helical rod. 

The pitch of the helix is the difference between the z-coordinate of the initial 
and final positions of one helical turn. A helical turn can be defined such that 
the projection of the vector position of the spatial curve along the xy-plane 
(vector xi of Eq. (13)), rotates of 2n around the z-axis. 

Fig. 3 shows g\(s) for the free standing helix of Fig. lb. We see that g\(s) 
oscillates, its maximum amplitude being smaller than 0.006. For the helical 
shape displayed in Fig. lc we found that the maximum value of g\(s) is smaller 
than 0.008 (data not shown). While for a circular helix g\ = 0, for the free 
standing helices displayed in Fig. lb and Fig. lc the function g x oscillates 
around zero with small amplitude. 

The small amplitude of these oscillations suggests that the relations TZ(s) ~ 
k F {s)[^k 2 F {s) + t 2 {s)]' 1 and V{s) ~ 2ttt f {s) (s) + r| (s)]" 1 , valid for cir- 
cular helices, could be used to derive approximate functions for the radius and 
the pitch of different types of helical structures, but the oscillatory behavior 
indicates that these relations are not simple functions of the geometric features 
of the helix. 
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6 Straight and planar inhomogeneous rods 

Straight rods (fcp = 0), and planar rods (/cp 7^ 0), have null torsion {jf = 0), 
and constitute particular cases of helices. In both cases there is at least one 
direction in space that makes a constant angle a = n/2 with the vector tangent 
to the rod centerline. 

The straight inhomogeneous rod is a solution of the static Kirchhoff equations 
that has non-constant twist density (Eq. (27)), in contrast with the homoge- 
neous case for which the twist density is constant. 

The twisted planar ring {k,F = Constant) is a solution of the static Kirchhoff 
equations only if the bending coefficient can be written in the form: 

B(s) = A cos(k F s) + B sin(k F s) + Cijk\ , (58) 

with Aq, B and Ci constant. If fcp is function of s (instead of being a constant) 
there exist no solutions for Eqs. (21-26). So, the existence of a planar solution 
related to the general form of the components of the twist vector given by 
equations (28-30) requires kp = Constant. 

7 Helical structure with intrinsic curvature 

The helical shape displayed in Fig. lb resembles that exhibited by the tendrils 
of some climbing plants. In these plants the younger parts have smaller cross 
section diameter, giving rise to non-constant bending coefficient. The main 
difference between the solution displayed in Fig. lb and the tendrils of climbing 
plants is that the solution in Fig. lb was obtained for an intrinsically straight 



21 



rod while the tendrils have intrinsic curvature [5]. 

The tendrils of climbing plants are stable structures while the helical solution 
displayed in Fig. lb is not stable because the rod is intrinsically straight [30]. 
We shall show that a rod with intrinsic curvature may have a static solution 
of the Kirchhoff equations similar to that displayed in Fig. lb. 

The intrinsic curvature of a rod is introduced in the Kirchhoff model through 
the components of the twist vector, k.(°\ in the unstressed configuration of the 



where k F (s) and Tp(s) are the curvature and torsion of the space curve that 
represents the axis of the rod in its unstressed configuration, simply called 
intrinsic curvature of the rod. We consider that the unstressed configuration 
of the axis of the rod forms a helical space curve with the intrinsic curvature 
satisfying 



where K and T are constant and B(s) is the bending coefficient of the rod. 
The linear constitutive relation (Eq. (18)) becomes 



M = B(s)(k 1 - fcHdx + B(s)(k 2 - k^)d 2 + C(s)(k 3 - k^)d 3 , (64) 



where C(s) is the twisting coefficient of the rod. The static Kirchhoff equations 
for this case, Eqs. (16), (17) and Eq. (64), are given by 



rod as 



4 0) =4 0) ( s )cos£, 

4 0) =£' + rf( S ), 



(59) 
(60) 
(61) 



B(s) k^(s)=K , 
B(s)rP(s) = T , 



(62) 
(63) 
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f[-f2h + f 3 k 2 


= 0, 


(65) 








= o , 


(66) 






fs -hh + h h 


= o , 


(67) 


(B(s)(ki 


-kP))'-B(s)(k 2 


-ki 0) )k 3 + C(s)(k 3 -ki 0) )-f 2 


= 0, 


(68) 


(B(s)(k 2 


- fc<°>))' + B(s)(h 


-k^h-Cis^h-k^ + f, 


= 0, 


(69) 




(C(s)(k 3 


-k^y + Bis^h-k^h) 


= . 


(70) 



The components of the twist vector are expressed as: 



k l = k F (s) sinx , 
k 2 = k F (s) cosx , 
h = x' + t f (s) . 



(71) 
(72) 
(73) 



In order to obtain the simplest solution for the static Kirchhoff equations Eqs. 
(65-70) with the intrinsic curvature given by Eqs. (59-61) and (62-63) we shall 
look for a solution such that x — £ m Eqs. (71-73). This solution preserves 
the intrinsic twist density of the helical structure. In this case, the Eq. (70) 
becomes simply [C(s)(t f — t^)]' = or 

M 3 = C(s)(t f - tP) = Constant , (74) 

and we obtain the following differential equations for the curvature k F (s), and 
the torsion t>(s), of the rod: 

[M 3 k F — B T F (k F - kf)\' - [B(k F - kf)]'r F = , 

[B(k F - 4 0) )]" + t>[M 3 k F — B T F (kp - 4 0) )] -f 3 k F = 0, ( 75 ) 

lB(kp-kP)]>kp + f> = 0, 

where we have omitted the dependence on s to simplify the notation. In order 
to obtain a helical solution of these equations we apply the Lancret's theorem, 
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Eq. (12), to the Eqs. (75). We obtain the following results: 



Ms) = , (76) 

[B(s)(k F (s) - kP(s))]> = =► k F (s) - kP(s) = JL , (77) 

[B(s)(r F (s) - 4%))]' = =► r F (s) - rP(s) = ^ , (78) 

where K and T are integration constants. From Eqs. (74) and Eq. (78) we 
obtain 

T = §gM 3 , (79) 



so that the ratio B(s)/C(s) has to be constant. 

From Eqs. (77), (78), (62), (63) and (12) we have 
k F (s) K + K 

— -= =tana. (80) 

Tp{s) T + Tq 

From Eqs. (68), (69), (77) and (78) it follows that 



A = 0, 
/ 2 = 0. 



(81) 



Therefore, the helical solution given by Eqs. (76-78) (obtained imposing x — 
f ) is a free standing helix (F = (f u f 2 , f 3 ) = 0). 

It follows from Eqs. (76-78), (62) and (63) that the solutions for the curvature 
k F (s), and the torsion t f (s), of the rod with helical intrinsic curvature are 
similar to those of intrinsically straight rods, Eqs. (41-43). Therefore, rods 
with intrinsic curvature and a non-constant bending coefficient given by Eq. 



24 



(48) (Eq. (49)) can have a three-dimensional configuration similar to that 
displayed in Fig. lb (Fig. lc). 

8 Conclusions 

The existence of helical configurations for a rod with non-constant stiffness has 
been investigated within the framework of the Kirchhoff rod model. Climbing 
and spiralling solutions of planar rods have been studied by Holmes et. al. [40]. 
Here, we have shown that helical spiralling three-dimensional structures are 
possible solutions of the static Kirchhoff equations for an inhomogeneous rod. 

From the static Kirchhoff equations, we derived the set of differential equations 
(31-33) for the curvature and the torsion of the centerline of a rod whose 
bending coefficient is a function of the arclength s. We have shown that the 
circular helix is the type of helical solution obtained when B(s) is constant, 
independently of the rod being homogeneous or inhomogeneous. 

Though the differential equations for the curvature and torsion are general, 
we have obtained only the simplest helical solutions (Eqs. (39) and (41-43)), 
obtained when the Lancret's theorem is applied to the differential equations. 
We show that these solutions are free standing and that the curvature and 
torsion depend directly on the form of variation of the bending coefficient. 
Figures lb and lc are examples of helical solutions of inhomogeneneous rods 
whose bending coefficients are given by Eqs. (48) and (49). The helical struc- 
ture displayed in Fig. lb is a conical helix since the projected curve onto the 
plane perpendicular to the helical axis is a logarithmic spiral, i. e., l/k\{s) is 
a linear function of s [9]. 
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In the particular case of an inhomogeneous rod with the intrinsic curvature 
defined by Eqs. (59-61) and (62-63), with B(s)/C(s) constant, we also ob- 
tain the helical solutions displayed in Figs, lb and lc. The tendrils of some 
climbing plants present a three-dimensional structure similar to that displayed 
in Fig. lb. In these plants, the cross-section diameter of the tendrils varies 
along them, giving rise to non-constant bending coefficient, and the differen- 
tial growth of the tendrils produces intrinsic curvature [5]. The bending and 
twisting coefficients of a continuous filament with circular cross-section are 
proportional to its moment of inertia /. It implies that B{s)/C{s) is constant 
for an inhomogeneous rod. Therefore, the tendrils of climbing plants can be 
well described by the Kirchhoff model for an inhomogeneous rod with a linear 
variation of the bending stiffness. 

Acknowledgements 

This work was partially supported by the Brazilian agencies FAPESP, CNPq 
and CAPES. The authors would like to thank Prof. Manfredo do Carmo for 
valuable informations about the Lancret's theorem. 

A Appendix: The differential equations for the curvature and tor- 
sion 

Here, we shall derive the Eqs. (31-33). Substitution of Eqs. (28-30) into Eqs. 
(21-26) gives: 

A - h (£' + t f ) + f 3 k FC os£ = , (A.l) 
+ fi (e + r F ) - h k F sine = , (A.2) 
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f 3 - h k F cosC + h k F sine = , (A.3) 

(S(s) /c_f sin 0' + (C(s) - B(s)) k F cosC (£' + r F ) - / 2 = , (A.4) 

(B(s) k F cosO' - (C(s) - B(s)) A; F sine (f + t f ) + /i = , (A.5) 

(C(s) {? + t f ))' = 0. (A.6) 

First, we extract f\ and f 2 from Eqs. (A.5) and (A.4), respectively: 

h = -(B(s) k F )' cose + [M 3 &f - B(s) k F r F ] sine , (A.7) 
f 2 = (B(s)k F )'sm£ + [M 3 k F - B(s)k F T F ] cost , (A.8) 

where M 3 = C(s) (£' + r F ) is the torsional moment of the rod that is con- 
stant by Eq. (A.6). Differentiating fi and / 2 with respect to s, substituting in 
Eqs. (A.f) and (A. 2), respectively, and using Eqs. (A.7) and (A.8), gives the 
following equations: 

{-(B(s) k F )" - r F [M 3 k F - B(s) k F r F ] + f 3 k F } cose + 

(A.9) 

{[M 3 k F - B(s) k F r F }' - r F (B(s) k F )'} sine = , 
{(B(s) k F )" + r F [M 3 k F - B(s) k F r F ] - f 3 k F } sine + 

(A.10) 

{[M 3 k F - B(s) k F t f }' - T F (B(s) kp)'} cose = . 

Multiplying Eq. (A.9) (Eq. (A. 10)) by sine (cose) and then adding the re- 
sulting equations, we obtain the Eq. (31) for the curvature and torsion: 

[M 3 k F - Bk F T F }' - (Bk F )'T F = . (A.ll) 

Multiplying Eq. (A.9) (Eq. (A. 10)) by - cose (+sine) and then adding the 
resulting equations, we obtain the Eq. (32): 

(B k F )" + k F r F (M 3 - Bt f ) - f 3 k F = . (A. 12) 
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Finally, the Eq. (33) is obtained by substituting Eqs. (A. 7) and (A. 8) in Eq. 
(A.3): 

(Bk F )' + & = 0. (A.13) 
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Fig. 1. Helical solutions of the Kirchhoff equations using the Lancret's Theorem, (a) 

circular helix solution for an inhomogeneous rod with constant bending coefficient 

B a = 1 (47); (b) and (c) Lancret helices for inhomogeneous rod with bending 

coefficient given by Eqs. (48) and (49), respectively. The parameters, in scaled units, 

are M 3 = 0.05, V = 0.9, and the total length of the rod is L = 130. k F {0) = 0.24 for 

the helical solutions displayed in panels (a) and (b), and kp(0) = 0.22 for panel (c). 

(a) (b) (c) 




Fig. 2. (a), (b) and (c) are projected curves of the helical solutions displayed in 
Fig. la, Fig. lb and Fig. lc, respectively. We used Eq. (13) to obtain the projected 
curves. 

0.006 



-0 .002 
-0 .006 




Fig. 3. g\(s) for the free standing helix solution displayed in Fig. lb. 
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